Quantum cellular automata quantum computing with endohedral fullerenes 
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We present a scheme to perform universal quantum computation using global addressing tech- 
niques as applied to a physical system of endohedrally doped fullerenes. The system consists of an 
ABAB linear array of Group V endohedrally doped fullerenes. Each molecule spin site consists of 
a nuclear spin coupled via a Hyperfine interaction to an electron spin. The electron spin of each 
molecule is in a quartet ground state S — 3/2. Neighboring molecular electron spins are coupled via 
a magnetic dipole interaction. We find that an all-electron construction of a quantum cellular au- 
tomata is frustrated due to the degeneracy of the electronic transitions. However, we can construct 
a quantum cellular automata quantum computing architecture using these molecules by encoding 
the quantum information on the nuclear spins while using the electron spins as a local bus. We 
deduce the NMR and ESR pulses required to execute the basic cellular automata operation and 
obtain a rough figure of merit for the the number of gate operations per decoherence time. We find 
that this figure of merit compares well with other physical quantum computer proposals. We argue 
that the proposed architecture meets well the first four DiVincenzo criteria and we outline various 
routes towards meeting the fifth criteria: qubit readout. 

PACS numbers: 03.67.-a, 03.67.Lx, 73.21.-b, 73.22.-f, 76.70.Dx 
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I. INTRODUCTION 



Quantum information processing (QIP), relies on the 
ability to perform a selection of unitary operations on 
a multi-partite system. More precise criteria that must 
be met by any physical implementation of a quantum 
information processor have been proposed in Q . In par- 
ticular, in a spin-based quantum computer architecture 
the spin sites (i) should be easy to physically manipu- 
late, (ii) easy to replicate in numbers, (iii) they should 
possess some type of inter-site interaction and (iv) they 
also should be somewhat isolated from their environment. 
Meeting such ideals in practice is very demanding. It is 
believed by some that a solid-state implementation could 
prove highly advantageous in terms of the scaling up of 
the number of spin sites and many solid-state propos- 
als for a quantum information processor are currently 
under study p, ffl H. However, to achieve (iv), usually 
requires the spin-sites to be located within a near-perfect 
crystal lattice and this leads to great difficulties in satis- 
fying (i) and (ii) above. If the interaction between spin 
sites is mediated by a predominantly electronic interac- 
tion, as is the case in ion-trap quantum computers and 
||, the quantum computer operation will be very sen- 
sitive to stray electric charges. Such systematic charge 
noise can be source of heating and may be very diffi- 
cult to eliminate. Additionally the use of local qubit ad- 
dressing is ubiquitous in very many quantum computer 
designs. However the use of such local qubit address- 
ing comes with a potentially very large overhead such as 
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the building of nanoscopic metal gates as in [Bj, or the 
greatly increased bandwidth required in NMR quantum 
computing in order to frequency address a large number 
of spins in a molecule Q. For small numbers of qubits, 
this overhead may be acceptable but the decoherence ef- 
fects of imperfect local addressing in a device containing 
many tens of qubits may be much more serious. Finally, 
the use of nuclear spins alone in a spin-based quantum 
computer design leads to a serious difficulty in achieving 
a sufficiently high initial spin polarization with which to 
initialise the quantum computer. This is a primary ob- 
stacle in liquid state NMR quantum computing and the 
problem remains in solid-state designs even when the nu- 
clei are cooled to milli-Kelvin temperatures. 

Faced with all these challenges, a number of us || |7|, |J, 
following ideas put forward by the author in p] , have con- 
sidered a molecular-based quantum computer which may 
offer a potentially more robust packaging of the qubit. 
The central idea is to encode a qubit into the spin sys- 
tem contained within an endohedrally doped fullerenc. 

In earlier publications we outlined the unusual 
properties of Group V C60 endohedral fullerenes and in- 
dicated their potential use in a solid-state implementa- 
tion. We, however, did not there give a detailed descrip- 
tion of a concrete quantum computing architecture for 
use with this material. 

In ||, a locally-addressed quantum computer design 
using endohedral fullerenes was proposed. Briefly this 
design stored the qubit on the nuclear spin of the endo- 
hedral atom and used the endohedral electron spins as a 
local bus. The interaction between neighboring qubits is 
mediated by the magnetic dipole force between neighbor- 
ing endohedral electrons. By constructing a linear chain 
of molecules and subjecting the chain to a large magnetic 
field gradient one can frequency address the electronic 
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spin of an individual molecule and through this perform 
universal quantum computation. 

This type of design has a number of serious drawbacks: 

(a) to achieve fast local addressing the resonant frequen- 
cies of neighboring endohedral electrons must be well 
separated. Thus the microwave bandwidth required to 
frequency address many frequency separated spins grows 
with the number of qubits and may become technically 
quite challenging to achieve beyond a few hundred MHz 
(which corresponds to a few tens of qubits) . In addition, 

(b) to be scalable one must perform quantum error cor- 
rection and to do so will require the simultaneous perfor- 
mance of various quantum gates within the device. Thus 
the number of independent spectrometer frequency chan- 
nels required to produce such multi-frequency microwave 
pulses will grow with the device. To engineer such a spec- 
trometer may be not be ultimately scalable. Finally, (c) 
to achieve local addressing in |^| requires the generation 
of a magnetic field which is highly stable in time and 
possesses a very large, and preferably very homogeneous, 
spatial gradient. Currently spectrometers are engineered 
to yield highly homogeneous, ultra-stable constant mag- 
netic fields (zero spatial gradients), and the engineering 
of magnetic fields desired in ||, on a microscopic scale, 
though possible, will require further development. 

The purpose of this work is to show that a quan- 
tum cellular automata quantum computing architecture 
(QCAQCA), can successfully be applied to operate with 
Group V C60 endohedrals. Further, the resulting design 
does not exhibit the above mentioned drawbacks seen 
in The QCAQCA design only requires an homo- 
geneous constant magnetic field. More importantly, in 
the QCAQCA, full scale universal quantum computation 
with quantum error correction only requires frequency 
addressing at a fixed, small number of frequencies, in- 
dependent of the size of the device. These two advan- 
tages mean that a quantum cellular automata based de- 
vice may be far easier to develop in hardware than any 
locally addressed architecture. 

For the design developed below we exhibit the spe- 
cific pulse primitives necessary for quantum cellular au- 
tomata quantum computation. We are also able to esti- 
mate rough figures of merit, or number of quantum gates 
that can be executed in a decoherence time. The esti- 
mated figure of merit obtained below compares well with 
other proposed quantum computer designs. 

In section II below we review the properties of X@G60 
{X = 15 A, 31 P), pertinent for quantum information 
processing. In section III we briefly outline the es- 
sential ingredients of QCAQCA. In section IV we de- 
velop a number of tools (pulse sequences), for manip- 
ulating the Hamiltonian of an alternating linear chain, 
ABABAB, of these endohedrals where A = 15 N@C60 and 
B = 31 P@C60. In section V we examine the possibility of 
an all-electronic QCAQCA and find that universal quan- 
tum computation seems not possible using the electrons 
alone. In section VI we consider an architecture which 



also includes the nuclear spins and show, by develop- 
ing new gate pulse sequences, that universal quantum 
computation can be achieved via a quantum cellular au- 
tomata design. We finally, in section VII, propose a num- 
ber of readout techniques that might be applicable to the 
proposed ensemble (a Type-II quantum computer design 
Q), and single-issue (or a Type-I quantum computer de- 
sign |J), quantum computer design. 



II. GROUP V ENDOHEDRALS 

The endohedrals 15 N@C60, 31 P@C60 and N@C70 all 
exhibit very sharp ESR spectra (l0| |llj]. This indicates 
the presence of free electrons within these molecules. It 
has been further shown, both theoretically and experi- 
mentally fl|, [if, 0, that the trapped atom (N or P), 
sits at the geometric center of the fullerene cage and 
that the electrons are in a S — 3/2 quartet ground 
state. The trapped atoms are extreme examples of com- 
pressed atoms []l2[ where the electronic wavefunc- 
tion of the trapped atom is repulsed away from the en- 
compassing carbon cage and suffers a spatial compres- 
sion. This unique type of "bonding" leads to significantly 
higher nuclear-electron wavefunction overlap and thus 
larger Hypcrfinc coupling constants than what is found in 
"free" atoms. For the case of C60 Group V endohedrals, 
the large symmetry renders this hyperfine coupling highly 
isotropic. Another consequence of this compression is 
the virtual lack of any electronic interaction between the 
trapped atom and the carbon cage. The trapped atom is 
motionally confined by a harmonic-like potential to the 
center of the cage. In almost all respects, the trapped 
atom behaves as a "free" (unbonded) atom, though spa- 
tially restricted to be within the fullerene cage. All of 
the above findings imply that such endohedrals behave 
as nanoscopic molecular neutral atom traps. 

The distributed ir bonding electrons on the C60 also 
act as an almost perfect Faraday cage, strongly isolating 
the electrons of the trapped atom from external electric 
fields [^6|. This latter observation implies that it would 
require the application intense local electric field gradi- 
ents ~ lV/nm, in order to alter the Hyperfine coupling 
constant of the trapped atom. Such an electrostatic ad- 
dressing scheme has been proposed to address nuclear 
spin qubits in the Phosphorus donors of the Kane design 
||, but with far lower field gradients. Such a scheme 
would prove very difficult to execute here due to the re- 
quirement of such intense electric field gradients. One 
very significant advantage of containing the spin site in- 
side a relatively large C60 molecule (diameter ~ lnm), 
is that such endohedrals can be nano-positioned using 
current scanning tunneling microscope techniques [fl7|| , 
and moreover neighboring C60 molecules can be nano- 
positioned on the Silicon(100)-2x surface to be as close 
as ~ I. lnm. The C60 molecules, on the Silicon(lOO)- 
2x surface, are fixed translationally at room temperature 
and also rotationally at lower temperatures. More specu- 
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latively, due to the electronic-wavefunction spatial com- 
pression experienced by these molecules, the chemistry 
of the doped Group V fullerene material is practically 
identical to undoped fullerene. It may then be possible 
to self-assemble large organised spin structures using ex- 
isting well-studied synthetic supra-chemistry techniques 

The very sharp ESR spectra from these molecules in- 
dicates very long relaxation times and more recent mea- 
surements have shown the electronic relaxation times are 
Ti ~ Is at T ~ 7°K, while T 2 ~ 20/zs, VT @. No 
nuclear relaxation times have yet been recorded but they 
are expected to be several orders of magnitude longer 
than the electronic relaxation times. The measured T 2 
time contains contributions from several sources, e.g. un- 
resolved dipolar couplings. The theoretical maximum 
of T 2 , in the complete absence of the unwanted inter- 
actions with other paramagnetic impurities, is the re- 
laxation time T\. Indeed, for the case of Phosphorus 
defects in isotopically ultra-pure 28 Si, the phase relax- 
ation time of the loosely bound Phosphorus electrons 
can be as long as T 2 > .1 msec with 1\ > 1 hour |l9| . 
There T 2 is limited by Hyperfine interactions with resid- 
ual 29 Si nuceli ^0|. The indicated increase of T 2 towards 
T\ in the situation of vanishing spin density is crucial 
for all current spin-based quantum computer proposals 
{§ |§ § |, % § Achieving such limits will be very 
challenging. Making these assumptions for the endohc- 
dral electrons we expect T 2 ~ T\ ~ lsec, at T ~ 7°K, 
with perhaps longer times at lower temperatures. Such 
Electronic T\ relaxation times (~ lsec), in such compli- 
cated molecules is highly unusual and may be unique in 
all non-crystal hosted spin-sites. Long relaxation times 
also add to the usefulness of such material as hosts for 
storing and manipulating quantum information. 

An essential ingredient for any quantum information 
processing is an inter-qubit interaction which can gen- 
erate entanglement. For Group V C60 endohedrals this 
cannot be a direct electronic exchange type interaction 
as the electronic wavefunctions are tightly compressed to 
be within the C60 cage. The C60's Faraday cage does 
not restrict magnetic interactions and neighboring endo- 
hedrals experience a significant magnetic dipole coupling, 
H D ~J(l-3 cos 2 9) [ZI Z S Z - 1- S] ~ J D I Z S Z , in the weak 
coupling limit. The strength of this dipolar coupling has 
been estimated to be Jo ~ 50MHz(l/r 3 ), where r is 
the separation between the neighboring trapped atoms 
in the endohedrals measured in nanometers |22| , and 9 is 
the angle between r and the external B field. The Hy- 
perfme coupling constants for 15 N@C60, 31 P@C60 have 
been measured [^3| , and are given in Table Q. 

From the above, to a first approximation, the full 
spin Hamiltonian for the pair AB, (A = 15 N@C60, 
B = 31 P@C60), can be given as 

H/h = g e [i e B z S A - g$n N Bj z A + A A S A I A 
+ g e » e B z S z B - g B N ii N B z lf + A B S B I B 
+ J D S A S B (1) 



where I(S), labels nuclear(electronic) spin, and we have 
made the secular approximation dropping terms which 
do not commute with the electronic Zeeman Hamilto- 
nian since uf armour = g e [i e B z > (J D , cj£ ormow = 
5jv/ijv-B z , A). The individual ESR and NMR spectra of 
these molecules are simulated in Figure [IJ and have been 
measured in [ pd[ . 

III. QUANTUM CELLULAR AUTOMATA 

Typical quantum computer architectures assume the 
capability of locally addressing every qubit. Such a 
requirement is highly challenging to engineer and the 
achievement of this level of individual quantum control 
is a major goal in almost every current implementation. 
The use of such local control methods forces an interac- 
tion of the quantum information with very many "classi- 
cal" control gates, each possibly providing a decohering 
effect on the quantum information. It is also important 
that the effects of such control gates on the computer's 
Hamiltonian remain as "classical" as possible in that they 
effect a change in the parameters appearing in the sys- 
tem Hamiltonian and do not themselves become entan- 
gled with the quantum computer. This criteria of classi- 
cal gating becomes non-trivial as the length scales of the 
quantum computer architecture reduces to nanometers. 

An alternative architecture is to utilize a small number 
of identifiable spins, placed in a regular spatial pattern, 
and to manipulate the quantum information encoded in 
this spin chain via global addressing techniques. This 
first of such globally addressed architectures was invented 
by Lloyd in 1993 p4j, and utilized three types of address- 
able spin arranged in the linear pattern ABC ABC ABC, 
and where each spin site encodes one logical qubit. In 
fact this cellular automata design was one of the very 
earliest quantum computer architectures proposed in the 
literature. Lloyd showed that such a quantum computer 
architecture was universal. Benjamin p5| , found a simi- 
lar architecture which used only two types of identifiable 
spin species arranged in the alternating linear pattern 
ABABAB. This reduction in spin resources came with 
an increase in logical encoding: a logical qubit is now en- 
coded into four spin sites with a buffer space of at least 
four empty spin sites between each logical qubit. 

The operation of Benjamin's architecture centers on 
the ability to implement the global unitary operator, Aj . 

Denoting the spin up(down) state as 1(0), A^ is the con- 
ditional application of the unitary U to the A spins in 
the alternating spin chain ABABAB, depending on the 
state of A's neighboring B spins. In particular, letting / 
be the sum of the states of the neighboring B spins, we 
have / € [0,1,2]. Thus, for example, AY is the condi- 
tioned application of U to all A spins in the alternating 
chain which have neighboring B spins that are different 
from each other, i.e. / = 1. One has a similar definition 
for for Bj , for the application of U on all B spins con- 
ditioned on their neighboring A spins. The case when 
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U ~ NOT, in AH (or Bj), occurs quite frequently in 
Benjamin's architecture and is shortened to Af (or Bf). 

Benjamin shows in [ p5| , that through a judicious se- 
quence of applications of A H and BH , one can implement 
single qubit operations and two-qubit CNOT operations 
and thus the architecture is universal. In particular, to 
move qubits through the spin chain one applies an al- 
ternating pulse sequence of A\ followed by B\ , while the 
generation of a control-U between two neighboring logical 
qubits requires ~ 30 global pulses ptjfl . One can trans- 
late all the standard circuit-based quantum algorithms 
to run with this quantum cellular automata architecture 
p6[ . This programming architecture, though somewhat 
expensive in terms of spatial and temporal resources, is 
ideal for use in systems where the individual control of 
qubits is difficult and it could provide a interim test-bed 
for various implementations while the more long term 
goal of local gating is developed. 

Since Lloyd's and Benjamin's pioneering work little has 
appeared in the literature with regard to quantum cellu- 
lar automata designs for quantum computers [p7| . In 
this paper we will consider a quantum cellular automata 
architecture consisting of an alternating linear array of 
Group V endohedrals, ABABAB, with the AB Hamilto- 
nian given by (Q) (schematically depicted in Figure ||). 
Essentially this task reduces to developing a "quantum 
algorithm" for generating Benjamin's global operation 
AH (similarly for BH). We first investigate the possi- 
bility of implementing AH , using the S (electron), spins 
alone. 



IV. THE GLOBAL OPERATION 

In this section we consider an alternating spin 1/2 
chain of electronic spins S A S B S A S B . We begin by con- 
sidering the S A S B unit alone where A spin has only one 
neighboring B spin. We assume that we know how to per- 
form the controlled operation of applying U to S A , with 
S B controlling, an operation we denote by C(B,A;U). 
We shall describe below the details regarding the se- 
quence of pulses required to execute C(B,A;U). We 
next consider the case when A has two neighbors, e.g. 
S B S A S B . Remembering that one cannot separately ad- 
dress an individual neighboring B spin one can show that 
instead of performing C(B,A;U), the above mentioned 
sequence of pulses now performs 

S(B, A; U) ee C(B, A; U)C{B 7 A; U) = C{B, A; U) 2 , 

(2) 

where Eq. (0), defines the function 3. One can observe 
that the case when the desired global operation is A\ = 
A^ OT , we have Ai = E(S,,A;NOT). 

We now indicate the procedure for implementing A 2 , 
when the desired U ~ exp(i9S A ). For ease of notation 
we define the standard NMR/ESR symbol for ir/2 pulses 
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[+5»]=exp(i|s») , (3) 

and the more general rotation through the angle 7 about 
the z— axis, for example, by the symbol 

{Z^}=exp(i 7 5^) . (4) 

We find that we can execute A 2 , when U ~ exp(i9S A ) 
through the following pulse sequence (read right to left), 

A u 2 = E(B,A;{X A })-{X A /2 _ s }-{Z A v/2 }-E(B,A;{X A }) 
■{Z A n/2 } ■ {X A 2 _ S } ■ E(B, A; {X A }) . (5) 

where 26 — 9. To implement A¥ one sets 6 = 9 and 
replaces the left most term in (0), E(B,A;{X A }), by 
{Z A }, while to implement An , one applies {X B } to all 
the B spins, then A 2 , and finally flip the B spins back 
with {X B ^}. The above pulse sequence (||), can be eas- 
ily modified to suit other forms of U besides exp(i9S A ), 
however we will find below that the most useful case is 
when = ir and U = exp(iirS A ) ~ NOT on the A spin. 

V. ALL ELECTRONIC QUANTUM CELLULAR 
AUTOMATA QUANTUM COMPUTER 

From Figure [j], we can separately address the elec- 
tronic spins S A and S B with a multifrequency selective 
soft pulse of length > 20ns. If we first consider both S A 
and S B to both be spin-1/2, then one can use standard 
pulse encodings of the CNOT |o), 

C(A,B,NOT) = [-S B ][TSfTS B }[±2SfSf}[+S B }e ±l ^\6) 
C(B, A; NOT) = [S A ] [ T S A T S B ] [±2S A S B ] [+S A ]e ± ^'{l) 

to generate the operation A\ and B\ of the previous sec- 
tion. To perform universal quantum computation one 
must be at least capable of implementing a two-qubit 
CNOT and a minimal set of one-qubit operations such as 
the phase gate, the 7r/8 gate, and the Hadamard gate. In 
the quantum cellular automata architecture of Benjamin 
this means that one must at least be able to implement 
the global operation A*j , where U is taken from this min- 
imal set of one qubit gates and the case where U ~ NOT. 
NMR and ESR pulses effectively generate terms in the 
Hamiltonian which arc linear in the spin operators, I x ,y,z, 
and S x ^ y>z . In particular, the spin operators when com- 
muted close to form the Pauli algebra. In the case of 
spin-1/2, the Pauli algebra generates the group SU(2), 
or all possible spin-1/2 unitaries. Thus one can gener- 
ate any desired spin-1/2 unitary through a sequence of 
pulses when acting on an isolated spin. Furthermore, in 
the case of the alternating chain ABABAB, one can find 
suitable alterations of (0), when used in (||) and (||), to 
generate any given unitary U in AH and similarly for K?. 
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Thus an all electronic quantum cellular automata quan- 
tum computer architecture (QCAQCA), using a spin-1/2 
ABABAB chain is possible. 

The endohedrals 15 N@C60 and 31 P@C60, however 
posses electron spin 5 = 3/2 and the application of Ben- 
jamin's QCAQCA is not as straightforward. One pos- 
sibility is to consider the four electronic levels of the 
electron spin with 5 = 3/2 to constitute a qu-dit with 
d = 4. Multi-level quantum computation has been pro- 
posed elsewhere |M, but to extend Benjamin's formula- 
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tion of the spin-1/2 QCAQCA to operate with qu-dits 
is far beyond the scope of this article. Indeed, as shown 
below there are reasons to believe that that such an ex- 
tension may not be possible within the confines of NMR 
andESR. 

Remaining with qubits we instead propose to encode 
a qubit into the four levels of a spin 3/2 in two ways: 
an "inner qubit" {|m s = ±1/2)}, and an "outer qubit" 
{|m s = ±3/2)}. By a small alteration of (^), one can 
obtain 



C inner (S A ,S B ;NOT) 



C outer (S A ,S B ;NOT) = {S®]{-S A - S®][+2Sf S®}[- S*}e+™ ' * 



C mner (S B ,S A ;NOT) = C outer (S Bl S A ;NOT) = [S A ][-S: 

I 



A 



S?][+2SfS?][-S?]e 



A|.±Wr/4 



(8) 
(9) 



where Ci nner is a CNOT on the Hilbert space spanned 
by the inner qubits, and similarly for C ou t er , S A = 
NOT(S A ), and we have used the spin-3/2 representa- 
tions of the SU(2) group. Armed with this one can easily 
construct the Benjamin global operations A\ and B\ on 
the inner qubit subspace. One can further construct (see 
section VI C| ), the global operations Aj OT , (and simi- 
larly Bj OT ) 1 for / = 0, 1 and / = 2, the operations 
most frequently used in the quantum cellular automata 
quantum computation. 

However one finds that one cannot generate all the 
one-qubit gates required for universal quantum compu- 
tation. As above, ESR pulses essentially generate terms 
in the Hamiltonian that are linear in the spin operators 
and again these terms commute to form a Pauli alge- 
bra. Evolution under such terms lead to dynamics that 
resides in the group generated by this Pauli algebra (or 
5/7(2)). However, since now we have 5 = 3/2, this 5/7 (2) 
is merely a subgroup of the full group of unitaries that 
can act on a spin-3/2, but as long as the system Hamilto- 
nian is linear in the spin operators we need only to con- 
sider the dynamics within the SU (2) subgroup. Normally 
a unitary operation U £ 5/7(2) C 5/7(4) will not factor 
into unitary operations on the inner and outer qubit sub- 
spaces. However if U is diagonal in the \m s ) basis it can 
so factor and with this one can find pulse sequences to 
implement Aj where U is a phase gate, and the ir/8 
gate but crucially not the Hadamard gate. Thus an all 
electronic-QCAQCA with a system Hamiltonian which 
is linear in the spin operators and where 5 > 1/2 seems 
not possible. 

One way of implementing 5/7(2) operations on the 
qubit subspaces and one which is well known in NMR 
and ESR is to introduce terms in the system Hamiltonian 
which are not linear in the spin operators. Such terms lie 
in the full spin-3/2 algebra and the generated unitaries no 
longer lie in the SU(2) subgroup but in the much larger 
5/7(4) group. One effect of such nonlinear terms, for ex- 



ample Zero-Field-Splitting terms: Hzfs oc 5f , can be to 
lift the degeneracy of the microwave transition frequen- 
cies between the \m s ) «-► \m s ± 1). This allows one to 
frequency address the transition \m a = —1/2) \m s = 
±1/2), or the inner qubit subspace, and directly imple- 
ment 5/7(2) operations on this subspace. This allows 
the generation of arbitrary one-qubit unitaries U [ |32| . 
However, this solution frustrates the coherence transfer 
portion, [+2S A Sf], of the pulse sequence (||), and (|]), 
and thus the construction of the global QCAQCA oper- 
ator Aj . This can be seen by realising that to obtain a 
[+2S A S®] pulse one must effectively cancel out all the 
terms in the system Hamiltonian bar the term S A S®. 
As standard in NMR and ESR one does this cancellation 
through "average Hamiltonian" theory p3| ], via the ap- 
plication of various 5/7(2) unitaries by RF or MW pulses. 
As long as the system evolution remains within the 5/7 (2) 
subspace such an averaging out can be achieved. How- 
ever when terms nonlinear in the spin operators are in- 
troduced into the system Hamiltonian no application of 
5/7(2) unitaries can, in general, average out the resulting 
"nonlinear" evolution (i.e. in the average Hamiltonian 
theory nonlinear terms in the Hamiltonian will almost 
always average to other nonlinear terms), and thus the 
generation of [+2S A S® ], and from this Aj , is not possi- 
ble. 



Thus we have found that when operating on an 
ABABAB, spin-3/2 chain with ESR pulses one cannot 
generate all of the necessary global operations required 
for universal quantum computation. We next find that 
the combination of nuclear and electronic spins within 
the chain can improve matters. 
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VI. NUCLEAR-ELECTRON QUANTUM 
CELLULAR AUTOMATA 

15 N@C60 and 31 P@C60 both have nuclear spin-1/2. In 
this section we assume that the quantum information is 
stored in the nuclear spins of the dopant atoms in the 
ABABA chain. This has the significant advantage that 
the nuclear spin relaxation times are typically longer than 
the electronic spins by several orders of magnitude. As 
the nuclei only are coupled via the Hyperfine interaction 
to the electrons, one can, via a suitable RF and MW pulse 
sequence, cancel out this Hyperfine interaction with high 
precision. This allows the nuclei to act as a quantum 
memory and store the quantum information for very long 
periods of time in between processing. 

To execute the Benjamin global operation AH we es- 
sentially use the inner electronic qubits as a local "bus" 
for the quantum information stored in the nuclei. The 
procedure can be thought of as a quantum algorithm for 
the generation of the global operation and consists of 
several steps which are pictorially represented in Figure 

& 

The algorithm starts with the inner electronic qubit 
bus in the ground state (| m s ) = | — 1/2) for all molecules), 
and this bus is returned to the ground state after im- 
plementing the global operation AH . We first (a) as- 
sume that we have an arbitrary pattern of quantum in- 
formation written onto the nuclear spins of the A and 
B molecules with the electrons set as above. We then, 

(b) swap the quantum information of all the B molecules 
from their nuclei to their inner electronic qubits. We then 

(c) tag those A molecules which will receive the unitary 
operation U in AH by flipping the state of the electronic 
qubit of all the A molecules conditioned on the state of 
it's neighboring B electronic qubits. We then (d) undo 
(b) by swapping the quantum information back into the 
B nuclei from their electrons. 

At this point all the inner electronic qubits are back in 
their ground states except for the "tagged" A molecules. 
We then, (e) perform a controlled-t/ operation on the 
nuclear qubits of all the molecules using their inner elec- 
tronic qubits as the control. Finally (f), we undo the 
operations (d), (c) and (b). The system is then ready for 
the execution of the next global operation. 

One restriction we make in the following is that we 
rule out the use of selective pulses on the nuclear spins 
of the Nitrogen and Phosphorus. From Figure |l|, to ex- 
ecute such selective pulses would require RF pulses of 
very long duration to frequency differentiate between the 
NMR transitions of the Nitrogen and Phosphorus. Such 
slow nuclear pulse sequences would lead to very long 
gates times. We will therefore not allow ourselves to use 
such selective nuclear pulse sequences and instead make 
use of fast "hard" nuclear pulse sequences Q. We do 
allow selective addressing of the electronic transitions of 
the Nitrogen and Phosphorus as these pulse sequences 
are much shorter in duration. With this restriction how- 
ever we are still able to design all the pulse sequences 



that are required to carry out the steps (b)-(e), above. 
The resulting combined total pulse sequence is somewhat 
lengthy but could be substantially shortened by applying 
the principles of optimal control theory |55) . 

We now expand on the above steps (a)-(e), in detail. 
In section (A) below we describe the initialization of the 
spin-chain. In section (B) we design the SWAP pulse se- 
quence required for step (b). In section (C) we design the 
electronic pulse sequence to implement the conditional 
operation required in step (c), and finally in section (D) 
we outline the execution of the controlled-?/ needed in 
step (d). In section (E) we discuss the undoing of the 
previous steps (b)-(d), while in section (F) we provide 
some estimates of the pulse durations and logical QCA 
gate times. 



A. Initialization 

The QCA architecture may be implemented on an en- 
semble (Type-II) or single-issue (Type-I), quantum com- 
puter. The achievement of high nuclear polarizations 
has proved to be a major obstacle in NMR quantum 
computation. It remains a source of difficulty here as 
well. Polarizing the electrons is far easier, and one can 
achieve a difference in the ground to excited electronic 
populations of e ~ 0.999, at a temperature of \°K and 
B ~ 10T. Even when multiple spins are tensored to pro- 
duce a pseudo-pure state, with such high individual po- 
larizations one can still achieve a pseudo-pure state with 
a purity on the order of e ~ .998, with 1000 spins at 
these temperatures and magnetic fields. This polariza- 
tion can be transferred to the nuclei via an INEPT pulse 
sequence |56| . To repump the electronic polarization and 
achieve a subset of completely spin (nuclear and elec- 
tron), polarized molecules one could consider two chains 
of equal length: an ABAB ■ ■ ■ chain and a CT>CT> ■ ■ ■ 
chain where all A, B, C, and V are globally address- 
able. Since at 1°K and B = 10T, half of the spins (the 
electrons), are completely polarised, one can consider an 
ABAB ■ ■ ■ ABABXCVCV ■ ■ ■ CVCV super-chain where X 
has a different resonant frequency yet again. One can 
then use a spin cooling quantum algorithm |}7j, to effi- 
ciently spin cool all the nuclear and electronic spins on 
one side of X (ie. X acts as a spin gate). One would then 
switch off the interaction between the two sub-chains at 
X and use the spin polarised half chain to perform the 
quantum computation. Such a design would be suitable 
for a Type-II (ensemble), quantum computer. Alterna- 
tively, in the case of a Type-I quantum computer one can 
replace the CVCV chain with an ABAB chain and place 
the readouts on the polarised half chain. More practically 
any Type-I quantum computer will require a single-spin 
readout of the electronic spins of the endohedral atoms. 
Such a readout may also serve to initialize the electronic 
qubits. We discuss possibilities for a single-spin readout 
further below. 
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Nuclear-Electronic SWAP 



A crucial ingredient in the above mentioned quantum 
algorithm to generate A^ is the Hyperfine SWAP. This 
operation performs a logical swap between a qubit stored 
in the nuclear spin and a qubit stored in the electron spin. 
Below we find two types of Hyperfine SWAP, one that 
performs the swap using the inner electronic qubit and 
one that uses the outer electronic qubit. For the most 
part we will use the former type of swap but we will find 
later that the latter swap, which incorporates the outer 
electronic qubit, may be of significant use in the problem 



J 



of qubit readout and we discuss this below. 

To develop the pulse sequence for the Hyperfine SWAP 
it is useful to break down the SWAP operation into three 
CNOT operations. Denoting SWAP (I : S), we have, 
SWAP(I : S) = CNOT(I;S)CNOT(S;I)CNOT(I;S), 
where CNOT(I; S) is the controlled NOT for I control- 
ling S. In our case 1 — 1/2 while S — 3/2. Substituting 
the spin-3/2 representations for S into the standard spin- 
1/2 CNOT pulse sequence @), fails to yield a proper 
CNOT operation. One instead must make reference to 
the two cases of a CNOT with respect to the inner and 
outer qubits of S. With this one can find: 



while 



With these one can derive 



CNOT outer (I;S) = [S y }[TS z Tlz}[±2S z I z }[-S y }e ±l ^ 4 , 

CNOT outer (S;I) = [I y }[±S z ±I z ][±2S z I z }[-I y }e ±l7l/4 , 

CNOT mner (I;S) = [S y ][TS z Tlz][±'2S z I z }[~S y }e ±t ^ A , 

CNOT mner (S;I) = {-I y ]{TS z Tlz}{±2S z I z }{I y }e ± ™/ 4 . 



-i-rr/4. 



SWAP mner (I:S) = [I v ] [+2I Z S Z ] [S y - I y ] [S x - 4] [+2I Z S Z ] [I x - S x ] [+2I Z S Z ] [I v ]t 
SWAP outer (7:5) = [-/„] [+2I Z S Z ][-I y - S y ] [-I x ~ S x ] [+2I Z S Z ] [S x + I x ] [+2I Z S Z ] [- I v ]e +3i ^ 4 



(10) 



(11) 



(12) 
(13) 



and the related pulse sequences are graphically shown 
in Figure |J. We note for later the curious case of (13), 
where we have in some sense "amplified" the magnetic 
signature of the qubit from a Amj = ±1, transition to a 
Ams = ±3, transition. In the following we will restrict 
ourselves to the inner electronic qubit subspace and thus 

©• 

Having found a SWAP operation between the nuclear 
and inner electronic qubits of either molecule we now 
must show how one can adapt the sequence (p"2|), to 
swap out only the B qubits while leaving the A nuclear 
qubits alone. We now consider the case of two coupled 
AB molecules with the Hamiltonian (Q). The desired 
selective swap action can be achieved by replacing the 
[+2I Z S Z ] terms in (|l2|), with the pulse sequence shown 
in Figure |. This new pulse sequence effectively averages 
out the Hyperfine coupling interaction in the A molecules 
thus turning off the action of the Hyperfine CNOT and 
SWAP gates for these molecules. The Hyperfine inter- 
action for the B molecules is left intact by the pulse 
sequence but is reduced in magnitude thus lengthening 
their Hyperfine gate durations. 

This particular pulse sequence has a number of advan- 
tages, the most important of which is that we do not 
selectively address the nuclear spins and thus can apply 
fast hard pulses to I A and I B simultaneously. We only 



take advantage of the selective addressing of the elec- 
tronic spins of A and B as outlined at the beginning of 
section [y| Alternatively, one can use soft multi- frequency 
nuclear pulses in the following to approximate the hard 
nuclear pulse. In more detail, from J38|, |39|, the action 
of this pulse sequence is to remove to all orders (within 
the secular approximation made in |n|), all interactions 
between S A ^ S B and S A «-> I A . This sequence also re- 
moves all the Zeeman terms in the AB Hamiltonian and 
one is thus left with H/h ~ A B S B I B . When inserted 
into (|l2l), the effect is to swap the qubit that was stored 
in I B , into the inner electronic qubit of S . As men- 
tioned above, the quantum information stored in I A is 
not swapped into S A , but these spins do receive the local 
unitary transformation, [S A ] [I A ] . The local unitary [S A ] 
is removed with a spin selective electronic pulse [— S A ], 
however we cannot remove the [1^] without performing 
a spin selective nuclear pulse, an operation we wish to 
avoid. We will see however, that this extraneous local 
operation on I A is not important for the tagging of the 
A molecules and furthermore the operation will be "un- 
done" instep (d) (c.f. Figured). 
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C. Electronic tagging of / 

Summarizing the previous steps: we have now swapped 
out the nuclear qubits of the B molecules in the 
ABABABAB chain into their inner electronic qubits. 
Additionally, all the inner electronic qubits of the A 
molecules are in their ground states \m s ) = | — 1/2), i.e. 
we are at step (b) in Figure [| We will now electronically 
"tag" those A molecules targeted by the global operation 
A^ . More precisely, we wish to flip the state of the A 
inner electronic qubit subject to the function /, with the 
control being A's nearest neighbor B's inner electronic 
qubits. We noted at the end of section [v|, that many 
inter-molecular electronic unitaries cannot be achieved 
with SU{2) pulses. However, the case of / = 1 is sim- 
ply the operation 2(S, A] NOT), (see section [V), where 
one uses (||) for the CNOT in the inner electronic qubit 
subspace. The case / = 2 is more difficult and one must 
make full use of the construction given in (||), to obtain 
the following pulse sequence 



S(fl, A; VWOT) ■ {Z A n/2 } ■ S(fl, A; VWOT) 
■{Z A , 2 }-E(B,A; ^fNOT)e 



3i7r/4 



(14) 



where E(B,A;VNOT) = C(S B , S A ; E) 2 , (see (§)), and 



C{S B ,S A ;E) = e 1 



k*S$ e -i\-K{S%+S'l-'lSiS%) e -i%-KS$ e i-K/4 

With this pulse sequence one achieves the required spin 
flip of the A molecule's inner electronic qubit up to a 
phase. This global phase is unimportant as it will be 
removed in step (e), when we undo the unitaries. Finally 
the case / = is dealt with by first flipping all B inner 
electronic qubits, then using the / = 2 operation above, 
and then flipping back. In the case of Figure [| we have 
chosen to implement the tag / = 1 operation in step (c). 
Once we have set the electronic tag on the appropriate A 
molecules we "undo" the SWAP, returning B's quantum 
information back into its nuclei and i?'s inner electronic 
qubits to the ground state, arriving at step (d) in Figure 



In (|l6|), f3, 7, and <5 are the Euler angles of the de- 
sired nuclear SU(2) operation, C{I,S;X) is the Hy- 
perfine control-NOT given in ([H]), the nuclear rotation 
operators are given by R$(9) = exp(— iO 5/2), where 
6 = x,y,z and standard pulse sequences for these spin ro- 
tations can be found in [p6[ . Further D = R z (/3)R y ('y/2), 
E = Ryi-yfflR^-iS + P)/2), F - R Z ((S - (3)/2), 
X ~ NOT, and a is a phase. Crucially one can show 
that DEF = I. From this one can see that the action 
of the sequence (|l^), on the nuclear spin is unity when 
the inner electronic qubit is not set. When the inner elec- 
tronic qubit is set however, the action of (|l^) is to execute 
U from ( p"7j ) on the nuclear spin up to the phase a. One 
then finally applies the local unitary exp(ia(S A + S?Y), 
to all molecules to yield the final phase factor a in (|1^). 
With this one has the very powerful capability of apply- 
ing any desired single-qubit unitary to the nuclear qubits 
of the tagged A molecules. However, the typical oper- 
ation in the QCAQCA is Af, and this can be executed 
more simply by a Hyperfine CNOT between the elec- 
tronic and nuclear qubits of all the molecules. In Figure 
||, we illustrate the application of A\ in subfigure (e). 



E. Undoing the unitaries 



Following the arguments made in section VI above, 

we execute step (e), that is: once we have implemented 

the desired unitary on the A nuclear qubits we "undo" 

steps (b), (c) and (d). We SWAP the nuclear qubit into 

the electrons on all B molecules, undo the tagging of the 

A molecules, and then SWAP the quantum information 

back into the nuclei of the B molecules. The system is 

left with all the inner electronic qubits in their ground 

states while the nuclear qubits have received the global 
\u 



Gate duration 



operation A 1 ! 



D. Implementing U 

We now can implement any desired unitary on the 
flagged A qubits. The construction of this operation is 
well known and consists of applying hard pulses 

p4j , to all of the A and B and nuclei interspersed by 
CNOT gates between the inner electronic qubits and nu- 
clear qubits of all the molecules. More precisely, the con- 
trolled execution of any given unitary operation on the 
nuclear spin / controlled by the inner electronic qubit S 
on all the molecules may be written as 



C(S,I;U) 
U 



DC(I,S;X)EC(I,S;X)F , (16) 
e la R z {P)Ry(l)Rz{5) 
= e ia DXEXF . (17) 



As we noted above, the resulting pulse sequence is 
quite lengthy. However it may be possible to compress 
many of the above operations. Also all nuclear pulses 
are fast hard pulses of little duration Q , while the pri- 
mary slowdown arises from the nuclear-electronic SWAP 
whose duration is limited by the value of the Hyperfine 
coupling constant of 15 N@C60 ( ~ 20MHz), and the sep- 
aration of the spectral lines of the A and B molecules 
in the ESR spectrum of Fig. @ (~ 50MHz). The typ- 
ical QCAQCA global operation, Af, entails a pulse se- 
quence in which the terms [+az^z], occurs 15 times, each 
with an average duration of 50ns, thus roughly bring- 
ing the cycle time of the global operation down to 1/xs. 
As noted above, the simplest quantum logic gate in the 
QCAQCA scheme requires ~ 30 global operations. This 
finally brings the cycle time for logical gates in the re- 
sulting QCAQCA to be roughly ~ 30[is. 
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We noted in section |J, on< 3 has T\ ~ Is at T — 7°K, 
and T 2 ~ 20/is for concentrated samples of 15 N@C60 
Group-V endohedral material. As mentioned above, re- 
cent experiments have indicated a dependence on concen- 
tration for T2 and it is commonly expected in all current 
spin-based quantum computer implementations that T2 
will rise towards the value of T\ in the limit of zero spin 
concentration. If T2 ~ Is then one can expect on average 
10, 000 — 30, 000, logical operations within this dephas- 
ing time (figure of merit). Indeed, at temperatures lower 
than 7°K one might be able to achieve far greater fig- 
ures of merit as the decoherence and dephasing times 
increase. A figure of merit of 10 compares very favor- 
ably with most alternative proposed implementations for 
quantum information processors. 



VII. READOUT 

One of the most challenging aspects of any solid-state 
implementation is that of readout. We first discuss the 
possibilities for an ensemble readout for the QCAQCA 
and later a single-spin readout for a Type-I implementa- 
tion. 

In the original scheme of Lloyd , the only spin sites 
which could be individually addressed were those at ei- 
ther end of the chain as these spins were frequency dif- 
ferentiated by having a single neighbor. This is also the 
situation here. However, if these two sites were the only 
readout sites on the chain then the architecture would not 
be scalable in the presence of noise. One ensemble read- 
out possibility would be to attach an electrically isolated 
paramagnetic adduct to every B molecule, ( Gd@C82 
may be a possibility). Such readout sites would typi- 
cally possess short relaxation times and one would have 
to isolate such sites from the operation of the primary 
processor using ESR pulses. It should be possible within 
the cellular automata architecture to selectively readout 
the state of a single qubit, that is, to transfer the state 
of a given logical qubit into the readout site. This can 
be done by adapting the single-qubit unitary operator 



U, pulse scheme in j25|. This pulse sequence is arranged 
so that a logical qubit (which is usually encoded in four 
spins, ABAB), gets driven into the state of a single spin 
which is then subject to a unitary operation, U, through 
the application of the global operation (say Brj , when the 
qubit has been forced into the spin state of a B molecule) . 
The spin neighborhood of this target B spin is arranged 
so that this global operation only effects the state of that 
spin alone with an appropriate choice for /. For read- 
out one could instead apply an operation B^ WAP , which 
swaps out the state of the spin into the readout site. This 
would be read and then SWAPed back. 

For a Type-I quantum information processor based on 
the ABAB chain of 15 N@C60 and 31 P@C60 one must 
be capable of performing a single electron spin readout. 
We suggest a number of possible technologies that may 
be capable of performing such a reado ut bel ow. Before 
mentioning these we return to section VI B and Equa- 
tion ( |l3| ) . There we noticed that this pulse sequence was 
able to SWAP the qubit between the highest and low- 
est weight subspaces of the nuclear and electronic spins, 
1=1/2 and S = 3/2, and in the process increase the 
detectable magnetic signature by a factor ~ 1000. 

We will now further enquire whether there exists sim- 
ilar pulse sequences which could further SWAP out the 
quantum information now stored in the outer electronic 
qubit into a separate coupled electronic system of larger 
spin. This could be very advantageous, for example, for 
the purposes of coupling into the spin-chain, a nano- 
molecular magnet with a spin of ~ 10 — 30 as a potential 
readout site Jig] . As the SWAP gate is built from CNOT 
gates we can focus on the latter. Further, the pulse se- 
quences must only consist of higher spin representations 
of the Pauli group. There may be many such pulse se- 
quences and we here present one set of sequences. The 
ones given below perform a CNOT between the highest 
and lowest weight subspaces of spins / and S. We set 
/ = a/2 and S = 6/2, where b > a and both a and 6 
are odd integers. Letting a = [a/2] and = [6/2], where 
[] is the integer part of the fraction, we can find CNOT 
pulses for the following four cases: 



a even, even : C\I,S;NOT) = [—S y ][—S z - I z ][+2S z I z ][S y ]e +i7T ^ 4 , (18) 

C ( S, I; NOT) = [-I v ] [-S z - I z ] [+2SJ Z ] [I y ]e + ™' 4 , (19) 

a odd, odd: C { I,S-NOT) = [-S y ][+S z + I z \[+2S z I z \[S y \e +m l A , (20) 

C { S, I; NOT) = [-/„] [+S Z + I z ] [+2SJ Z ] [+I v \e+"' i , (21) 

a odd, even: C { I, S; NOT) = [S y }[+S z + I z \[+2S Z I Z ][-S y ]e+"' 4 , (22) 

C ( S, I; NOT) = [I y ] [-S z - I z ] [+2SJ Z ] [~I y ]e+™l 4 , (23) 

a even, odd : C\I, S; NOT) = [S y ][-S z - I z ][+2S Z I Z ][-S y ]e + ™/ 4 , (24) 

C { S, J; NOT) = [I v ][+S z + I z ][+2S z I z ][-I y }e +t ^ 4 . (25) 

(26) 
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A Type-I quantum computer design will require the ca- 
pability of single qubit readout and this translates in our 
case to the capability of reading out the spin state of the 
cndohedral electrons. One can classify the various read- 
out methods into the general categories (i) force, (ii) elec- 
tric, and (iii) optical, measurement techniques. As was 
mentioned before, since the electronic wavefunction of 
the trapped atom is totally confined within the C60 cage, 
spin measurement techniques that involve the physical 
transport of this spin outside the molecule are not pos- 
sible here. Some techniques, such as ODMR (optically 
detected magnetic resonance) , have so far not proved pos- 
sible with 15 N@C60 and 31 P@C60, and thus the addi- 
tion of separate readout sites which can be dynamically 
coupled and decoupled from the primary processor via 
NMR/ESR pulses are warranted. 

Possible techniques in the above mentioned cate- 
gories include (i) Magnetic Resonance Force Microscopy 
(MRFM) HI]; (iia) Micro-Squids, a technology which is 
already capable of discriminating a Ams = 30 |44| ; (iib) 
Scanning- Tunneling-Microscope ESR, a technique which, 
though not well understood, has yielded single molecule 
ESR spectroscopy of iron impurities in Silicon , and 
in single BDPA molecules on Silicon [Q. (iic) A sin- 
gle molecule endohedral fullerene single electron transis- 
tor p7| . This is a very challenging technology which 
has yielded a single electron tunneling current through 
a C60 molecule with an electro-mechanical coupling to 
the quantized motion of the entire molecule. Performing 
such an experiment with 15 N@C60 may yield endohedral 
spin information; (iiia) Magnetic coupling of the endohe- 
dral electronic spin to a solid-state paramagnetic optical 
dipole such as a nanocrystaline N-V center in diamond. 
The optical paramagnetic center can then be probed 
via optical shelving techniques pq ]; (iiib) ODMR via a 
magnetic coupling to a paramagnetic endohedral adduct 
which is chemically bonded to the 15 N@C60 molecule. 
For use with ODMR the adduct endohedral should ide- 
ally possesses an optical transition in the visible spectrum 
such as Er^Scg-xNiaCSO @ . 

The techniques suggested above are, in some cases, al- 
ready at the single electron detection level (though not 
capable of detecting the spin orientation), while the out- 
coupling of the quantum information into a large-spin 
system may allow the techniques of micro- Squids and 
MRFM as they presently stand to act as a readout of 
a single endohedral electronic spin. 



VIII. CONCLUSION 

The philosophy taken in this work with regard to the 
construction of a quantum information processor is to fol- 
low as closely as possible the system Hamiltonian that na- 
ture provides. The globally addressed approach has very 
significant advantages over the local gating approach. 
Global addressing avoids the buildup of numerous deco- 
herence pathways associated with the effects of the local 



gates on the processor. Global addressing also avoids the 
very significant problem of the scaling up of the external 
resources required to execute local gating, eg. numerous 
metallic contacts or multi-frequency MW/RF generators 
with increasing bandwidths. Although it may be true 
that for local gating, the external resources required in- 
crease polynomially with the number of qubits in the 
processor, such an overhead may not be experimentally 
feasible as the size of the processor grows large. 

In some sense the technical complexities involved in 
the construction and operation of local gating are trans- 
ferred into the quantum cellular automata "software" . As 
shown above, the detailed operation of a quantum cellu- 
lar automata quantum computer architecture does not 
follow the standard quantum circuit model. This should 
not be seen as a disadvantage as all quantum circuit algo- 
rithms can be "compiled" to run on a QCAQCA with a 
polynomial overhead f^j . Indeed the QCAQCA may be 
able to run programs that do not follow the quantum cir- 
cuit model (see 50 for an example of a non-circuit quan- 
tum program in another architecture). Furthermore, as 
the fundamental quantum operations are implemented 
via the very well developed methods of NMR and ESR, 
the fidelity of software execution is heightened. The 
QCAQCA also has the very important property that it 
can execute quantum computations in parallel, a charac- 
teristic necessary for the ultimate scalability of the archi- 
tecture. 

In this work we considered the endohedral fullerene 
material 15 N@C60 and 31 P@C60 and argued that the 
system Hamiltonians and interactions present are suffi- 
cient to implement universal quantum computation via 
a quantum cellular automata architecture. Following the 
DiVincenzo criteria: (i) we argued that the Group-V en- 
dohedral materials behave as essentially electrostatically 
isolated nanoscopic molecular neutral atom traps that 
can be physically manipulated with relative ease using 
STM techniques and perhaps may be self-assembled us- 
ing fullerene supra-chemistry. (ii) the electronic decoher- 
ence times are perhaps the longest seen in any molec- 
ular system, of the order of seconds at temperatures of 
~ 7°K, which is also an upper limit for the electronic 
dephasing time. The nuclear decoherence times have not 
yet been measured but they are expected to be several 
orders of magnitude longer than the electronic decoher- 
ence times, (iii) The complete polarization of both the 
nuclear and electronic computational spins is feasible as 
half of the spins (electrons), within the entire system are 
completely polarised at moderate conditions of tempera- 
ture and magnetic field strength. There are efficient cool- 
ing schemes exist to shuffle the unpolarized spins away 
from the computational spins, (iv) entanglement can be 
generated via the inter-molecule magnetic dipole inter- 
action, the strength of which has been measured to be 
~ 50MHz. Armed with the two species 15 N@C60 and 
31 P@C60, we showed that one has sufficient control to 
implement a two-component quantum cellular automata 
quantum computing architecture storing the quantum in- 
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formation in the nuclear spins while using the electrons 
as a local bus. We further found that the expected "fig- 
ure of merit" compares very well with other proposed 
solid-state quantum computer designs. We also discov- 
ered pulse sequences that SWAPed the quantum infor- 
mation between the highest and lowest weight spaces of 
two spins of different size, e.g. S = 1/2 «-> S = 11/2. 
These sequences could prove useful for out-coupling the 
quantum information into a spin readout system with 
large spin. 

DiVincenzo's final criterion: (v) an efficient single 
qubit readout, is perhaps the most difficult challenge 
for any solid-state based quantum computer design. We 
have proposed various possible ensemble and single-issue 
(Type-I), readout technologies many of which are them- 
selves the subjects of intense study. The ideas pre- 



sented here combine the tremendous power of NMR and 
ESR science together with the very "clean" , and almost 
atomic, systems presented by the endohedrals 15 N@C60 
and 31 P@C60. These, married with supra-molecular 
chemistry, may provide a very real possibility for a phys- 
ical implementation of a quantum information processor. 
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B z = 2T 




A = 15 N@C60 


B = 31 P@C60 


Electronic Zeeman Energy 


g e ^eB z 


56GHz 


56GHz 


Nuclear Zeeman Energy 


gN/^NBz 


-6.1MHz 


34.5MHz 


Hyperfine Coupling Constant 


A 


21.2MHz 


138.4MHz 



TABLE I: Table of Hyperfine coupling constants and interac- 
tion energies for interacting 15 N@C60, 31 P@C60 molecules at 
B z = 2T. 
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FIG. 1: Simulated NMR (in subfigure A), and ESR (in sub- 
figure B), spectra of A = 15 N@C60 and B = 31 P@C60 with (a) 
and (b) labeling the spectral lines from A and B respectively. 




FIG. 2: Schematic depiction of the ABABAB quantum cellu- 
lar automata chain using the Group-V endohedral fullerenes. 



B A B A B A 




FIG. 3: Schematic representation of the steps (a)-(e) de- 
scribed in section |^j to execute the global operation A'f . 
The alternating molecular chain BABABA is depicted with 
the inner electronic or S qubits (larger shapes) , connected to- 
gether via the magnetic dipole interaction, while the nuclear 
or / qubits (smaller shapes), are connected to the electronic 
qubits via the Hyperfine interaction. The subfigures depict: 
(a) The S qubits in their ground state with some pattern of 
/ qubits; (b) A Hyperfine SWAP of the quantum informa- 
tion in the B molecules; (c) Conditional flipping of the A 
electronic qubit depending on the neighboring B electronic 
qubits; (d) Hyperfine SWAP of the B qubits back into the 
nuclei; (e) Controlled-(7 applied to the nuclei of all molecules 
conditioned by the state of their electronic qubits. Here the 
effect of AH is to flip the state of the leftmost A nuclear qubit. 





FIG. 4: Pulse sequence for a logical SWAP between spin I = 
1/2 and the inner (subfigure (A)), and outer qubits (subfigure 
(B)), of spin S = 3/2 as described in @ and @. The solid 
black bars are tt/2 pulses while the white bars are refocusing 
7r pulses. 
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FIG. 5: Pulse sequence to selectively couple / with 5* while 
decoupling S A with S B , and simultaneously decoupling S A 
with I A . The black pulses are 7r pulses and the duration 
of the entire sequence is matched to yield [+2lf Sf ] in (|l2[). 
The gray pulses are fast hard pulses which are required 
to extend the duration of the [2/ z 5 z ] for the case when swap- 
ping between the nuclear and electronic spin in 31 P@C60. 
By placing these pulses at the points (m — l)/(4m), and 
(3m + l)/(4m), in each third section as shown, we can ex- 
tend the pulse sequence by a factor m. Setting m = 6 brings 
the 31 P@C60 SWAP to be ~ 20MHz, and thus long enough 
for the electronic spin selective pulses to operate on S A and 
S B . 



